Microlensing of Kepler Stars as a Method of Detecting Primordial Black Hole Dark 
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If the Dark Matter consists of primordial black holes (PBHs), we show that gravitational lensing 
of stars being monitored by NASA's Kepler search for extra-solar planets can cause significant 
numbers of detectable microlensing events. A search through the roughly 150,000 lightcurves would 
result in large numbers of detectable events for PBHs in the mass range 5 x Mq to 10"'* Mq. 

Non-detection of these events would close almost two orders of magnitude of the mass window for 
PBH dark matter. The microlensing rate is higher than previously noticed due to a combination of 
the exceptional photometric precision of the Kepler mission and the increase in cross section due to 
the large angular sizes of the relatively nearby Kepler field stars. We also present a new formalism 
for calculating optical depth and microlensing rates in the presence of large finite-source effects. 

PACS numbers: 95.75.De, 95.35. +x,98.35.Jk 
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PRIMORDIAL BLACK HOLE DARK MATTER 



Primordial Black Holes (PBHs) have been considered 
as a candidate for dark matter (DM) since the days of 
Hawking [l|, Q and PBHs are recently re-emerging as ob- 
jects of intense study (see references [1, ^3] for recent re- 
views of constraints). PBH DM can form during the 
early universe as a result of large density fluctuations 
coming from, for example, bubble collisions, collapse of 
string loops, features in inflation potentials, etc. (see 
reference[5| for more examples and references). If PBHs 
form early enough they can evade the Big Bang Nucle- 
osynthesis limits on baryons, satisfy CMB constraints, 
and make up the entirety of the dark matter. 



While there is no compelling mass or mass range for 
PBH DM, there has been extensive experimental and 
theoretical work that has eliminated most mass ranges, 
starting with Hawking's famous limit, mpBH > lO^'^^g 
(mpBH > 1Q~^^Mq) for PBHs to not have evaporated 
by today. See reference Q for a recent summary of these 
constraints. Currently most mass ranges from 10~^*Mo 
to IO^^Mq are ruled out, with the exception being the 
flve orders of magnitude between 1Q~^^Mq < mpBH < 
which we call the PBH DM window. As we 
show below, microlensing of Kepler source stars can de- 
tect or rule out PBH DM over a large fraction of the 
PBH DM window. Thus, if the DM consists of PBHs, 
the experiment we propose here has an excellent chance 
of detecting it. Since gravitational microlensing is sen- 
sitive to any massive compact halo object (MACHO) a 
detection or limit would actually be more general than 
just PBH DM. 



KEPLER MISSION 

The NASA Kepler satellite is a 1 m aperture telescope 
with a 115 deg^ field-of-view in an Earth trailing helio- 
centric orbit, which is currently taking photometric mea- 
surements of around 150,000 stars every 30 minutes 6-8]. 
The telescope was launched in March 2009 and will point 
at the same fleld (in the Cygnus-Lyra region of the sky) 
for at least 3.5 years. The goal of the Kepler mission is to 
flnd extra-solar planets via small decreases in stellar flux 
due to rare planetary transits in front of the host star in 
edge-on systems. Very precise photometry is required to 
detect the tiny decrease in flux that an Earth size planet 
causes during a transit. 

While much of the data is still proprietary, in what 
follows we analyze a small portion of it to investigate 
its potential use in a microlensing experiment JeJ] . Look- 
ing through a portion of this publicly available data we 
see that Kepler source stars have V magnitudes between 
roughly 10 and 16, have distances from the Earth be- 
tween 0.9 and 3 kpc (with most between 0.9 and 1.6 
kpc), have stellar radii, i?,, between 0.9i?o and l.Si?© 
(for main sequence stars, but there are substantial num- 
bers of giant stars as well with 5Rq < i?, < 2Oi?0), and 
have photometric errors per observation between 20ppm 
and lOOOppm (with most between 300pm and 1000pm). 



MICROLENSING AND ANALYTIC ESTIMATE 

Microlensing searches for low mass objects have been 
performed and have returned signiflcant limits on any 
massive compact halo objects, including PBHs, that 
might constitute the dark matter (lol - [l3| . 

These microlensing surveys examined tens of millions 
stars in the Large Magellanic Cloud (LMC), which is at a 
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distance of 50 kpc, for periods of many years. Since the 
naive microlensing optical depth increases with source 
star distance and the total microlensing event rate is 
proportional to the number of monitored stars times the 
duration of the survey, it might be surprising that the 
many fewer Kepler lightcurves on much closer stars, can 
provide stronger limits on low mass DM in certain mass 
ranges than these surveys have. 

As calculated below, this is due to several factors. 
First, the extreme precision of the Kepler photometry 
allows very small magnifications to be detected. Re- 
quiring a sequence of 2-sigma or 3-sigma measurements 
on a typical Kepler lightcurve allows a magnification 
threshold of At = 1.001 or lower to be set (magni- 
fication, A = (Flux/ (Flux)) — 1). Since detection is 
said to occur whenever the magnification A > At, the 
effective microlensing tube[3] is wider by a factor of 
UT = \2At{A\ ~ l)-i/2 _ 2]i/2 ^ 6. 



Here we are using notation[14l Il6l | in which the 
Einstein ring radius (radius of image ring arising 
from perfectly aligned source and lens) is r^; = 
0.0193v/2:(l - x)[(i:'s/kpc)m9]i/2i?Q, where mg = 
(to/IO^^Mq), is the mass of the PBH DM lens, Ds is 
the source star distance, and x = Di/Dg is the ratio 
of the lens distance to source distance. We define the 
standard microlensing "tube" as that volume along the 
line-of-sight inside the Einstein ring, and define u to be 
the transverse distance between the lens and the source 
center in units of r^. The important parameter, ut, is 
the "effective" radius of the microlensing tube in units 
of rE in the point source limit; thus the usual standard 
mt = 1 implies At = 1.34. 

The usual point source (PS) optical depth (total 
number of PBHs inside the microlensing tube), Tpg — 

dxDsp{x)nuT'''^/'ni, is thus larger than in the stan- 
dard case by a factor of u^. where p{x) is the run of DM 
density along the line-of-sight (los). 

However, just as important for low mass PBH mi- 
crolensing is the fact that the effective microlensing tube 
radius is not determined by the Einstein ring radius, but 
by the radius of the source star due to finite-source ef- 
fects. For very low mass PBHs, the Einstein ring ra- 
dius is much smaller than the projected source star ra- 
dius, xR^,, and ut in the formula for optical depth should 
be replaced by uthicsh which is the value of the source- 
lens transverse separation that gives A — At- The 
quantity Uthicsh niay be approximated by = xRi,/rE, 
though below we make a better approximation using the 
finite-source lightcurve formulas from Witt & Mao [15]. 
For m « 10~^Mq with x « 0.5 and a typical Kepler 
source star with R,, ~ ^Rq and Ds ~ 1 kpc, we find 
M* = 51.9, which gives an optical depth (and total event 
rate) around 2,700 times larger than naively expected. 
(See below for a more accurate estimate.) Similarly, the 
usual statement that the optical depth is independent 
of DM mass is not true in this case. Naively, a lower 



mass means a narrower microlensing tube, but more DM 
objects, effects that cancel each other out in the opti- 
cal depth. For the case where all events are finite-source 
events, lowering the mass does not change the tube ra- 
dius, so the optical depth is roughly inversely propor- 
tional to mass, increasing the sensitivity to much lower 
mass DM objects. Likewise, the average event duration, 
(^e), defined as the time for which the source star is mag- 
nified above At-, is now roughly independent of mass, and 
the total event rate, F — t / {t^), now increases roughly 
inversely with decreasing mass. 

Of course, finite-source microlensing has been studied 
since the paper by Witt & Mao[l^, who show that in 
the finite source regime the magnification is limited to be 
less than A,„ax = yl+l/uf. For low mass PBHs, the 
lightcurve jumps quickly from near A = 1 to A = Amax 
and then stays roughly constant (actually following the 
limb-darkening shape of the source star) as the small Ein- 
stein ring transits the stellar limb, jumping back to near 
A — 1 ai the edge. See references [lOl, ll5| for example 
lightcurves, which are quite different from standard mi- 
crolensing lightcurves. The very precise Kepler photom- 
etry is what allows these low magnification finite-source 
lightcurves to be detected. For a given At, we use the 
previous equation to find a cutoff, M*max, such that a 
microlensing event is observable whenever u is less than 

'W^max — ^1 \j -^T 

The total optical depth and the total event rate in- 
clude integrals along the los, but the limits to these in- 
tegrals need modification when considering strong finite- 
source effects. As a; — J- 0, both the projected source star 
radius, xi?*, and te, go to zero. Their ratio, which 
is u», also goes to zero as x — 0, while as a; 1, 
rE ^ and ^ co. Because of this last depen- 
dence, Aniax — 1 as X — 1, and so there is always 
some limiting value of distance, Xmax, where only ob- 
jects that enter the tube with x < Xmax give detectable 
magnifications. We find Xmax ~ 1/(1 + where 
K = 51.9(ii»/i?0)/(u*maxV^"^^57/kpc). Thus the op- 
tical depth integration limits above should be between 
X — and x — Xmax, not between and 1. 

For the pure finite source case (u* 3> mt) and a con- 
stant DM density, the above optical depth integral can 
be easily performed. We find rps ~ spR^Xi^^x^ / m k, 
4.2 X 10-6a;„,ax^(i?,/i?o)'(^s/kpc)/m9. For At - 
1.001, i?* = li?0, Ds — 1 kpc, and mg — 1, we have 
a^max = 0.426 and rps = 3.2 x 10"''. To estimate the 
total event rate in this case we can use the very approx- 
imate formula from reference [lij or 14 1, replacing utte 

with Xmax-R*- We find Fpac-FS ~ 2TFSt'c/(7rXinax^*) = 

4.8 X 10^'^ events/year/star. The corresponding average 
event duration is given by {te) — t/F « 0.59 hours. This 
is a remarkably large event rate and if 150,000 stars were 
followed it would result in a total of around 720 events 
per year. 

However, another complication is that for a microlens- 
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ing event to be detected, it must contain enough signifi- 
cant measurements to be unlikely to occur by chance. For 
example, given Kepler's 30 minute cadence, monitoring 
150,000 stars means about 2.6 billion flux measurements 
per year. Assuming Gaussian errors, the probability of 
finding 4 sequential measurements 3-sigma above average 
is such that less than one such instance will occur in 3.5 
years of data. This requirement implies an event dura- 
tion, te, the time for which A > At, of at least 2.0 hours. 
This means that we cannot use the total event rate F, 
but must integrate the differential event rate ^ over a 
relevant range of event durations. This differential event 
rate is given in Equation 17 of reference 0| for non- finite- 
source events. In the limit of large this formula can be 
used just by replacing ut with m* throughout. There are 
also two simplifications to this formula due to the unique 
position of the Kepler field. First, the Kepler field is 
at galactic longitude and latitude {l,b) = (76.32°, 13.5°) 
which places it just out of the plane of the Milky Way 
disk, in almost the same direction in which the solar sys- 
tem is moving The Kepler source stars, at distances of 
1 to 3 kpc, are thus at nearly the same distance from the 
Milky Way (MW) center as the Sun (8.5 kpc). Thus, un- 
like the case for the LMC, there is no need to model the 
halo density and one can just use the local dark matter 
density p « 0.3 GeVcm"^ = 7.9 x IO'^Mqpc"^ along the 
entire los. Second, since the Kepler stars are in the di- 
rection of solar motion and also orbiting the MW center, 
there is no additional large transverse velocity of the Sun 
or source to include. Thus Equation 17 of reference 14 
becomes 



dr 



P n 2 



dxl3'giP'), 



(1) 



where g{p') = dyy'^^^il - y)-^/^e-^'y = 

§e-^'/^[Io{/3'/2) - (1 + l//3')/i(/372)], lo and h 
are modified Bessel functions of the first kind, 
/3' = 4rK_,/(t2„2)^ y = vl/{(3'vl), V, « 220 
km/s, is the halo circular velocity, and we introduced 
Wthrosh which stauds for ut in the case of point source 
lensing, in the case of pure finite-source lensing, and 
below will be approximated as something in between for 
the general case. For use later we note that ^(O) — |7r, 
and g{^') - 3^/?'-5/2 ^^^^^ ^^^j^^g ^, ^ 



NUMERICAL ESTIMATE 

In order to make a reasonably accurate estimate of 
the potential sensitivity of a search through the Kepler 
lightcurve data we need to integrate Eq. [1] from = 
'train to somc reasonable upper limit. We also want to 
use realistic distributions of Kepler star distances, radii, 
and At, and we want an approximation for Uthresh that 
interpolates accurately between Wthresh ~ mt at x near 



0, and Mthresh ~ M* at o; near 1. Ideally, we would also 
include the effects of limb- darkening, but do not do this 
here. We have done preliminary exploration of this effect 
by creating limb-darkened finite-source lightcurves and 
calculating wthrosh for these, but do not see a large effect 
for the few PBH mass values we checked. 

We looked at a subsample of around 5000 publicly 
available Kepler lightcurves from the NASA MAST 
websiteQ. We used data from the third publicly avail- 
able quarter and, besides the lightcurves of fluxes and 
flux errors, we have for each star: the stellar radius, i?*, 
Sloan r and g magnitudes, effective temperature, Teff, 
star position, extinction parameters Ay and E{B — V), 
etc. We estimate the apparent visual magnitude, V — 
g - 0.0026 - 0.533(5 ^ r) and the stellar distance 
from Ds = 1.19xlO-3i?,(reff/T0)2loO-2(^-'4v'+B.C.)]^pc^ 

where B.C. is the bolometric correction. Note that we 
make a very crude bolometric correction, using only the 
effective temperature and whether the source is a main 
sequence or giant star[l8|, but we include it because it 
slightly reduces the distances to the sources, thereby re- 
ducing the expected detection rate, and we want our cal- 
culation to be conservative. 

To find At for each source star, we calculated the aver- 
age of the reported flux errors over 300 data points near 
the middle of each lightcurve, and then estimated each 
At as one plus three times this average (for a 3-sigma 
detection requirement). We also calculated the standard 
deviation of the flux measurements over this portion of 
the lightcurve and found reasonable agreement with the 
average flux error. For our subsample of 5000 stars we 
find 2 X 10~^ < At < 3 x 10"'^, with vast majority having 
At ~ 0.001. 

As noted above, the approximation, Wthresh ~ 
misses some events, since for nearby lenses the projected 
source radius becomes smaller than the Einstein ring ra- 
dius, and therefore itthresh should approach ut rather 
than zero. To develop a better approximation, we numer- 
ically calculated a large set of finite-source lightcurves us- 
ing the formula from Witt and Mao[l5| , and then for each 
lightcurve calculated the actual value of Uthrosh (value 
of lens-source transverse separation for which A > At) 
as a function of u^,/ut- We find that Uthresh/ur is a 
fairly universal function of u^/ut, and therefore we can 
calculate for each value of x and find uthicsh using 
this universal function. Our fit function is: uthrcsh ~ 
mt(1 + ■47(u*/mt)^), for u*/ut < 0.75, Uthrcsh ~ 
for u^/ut > 4.5, and uthrcsh ~ w*/(Z]i Ci(M*/MT)*), with 
Co = 0.0971, ci = 0.925, C2 = -0.384, cg = 0.0723, and 
C4 — —0.0051, for 0.75 < u^/ut < 4.5. This approximate 
Mthrcsh differs from our numerically calculated values by 
no more than 2.5% over the mass range we investigated. 

We then numerically performed the 2-dimensional inte- 
gral of Eq.[T]for various values of PBH mass, various S/N 
requirements (e.g. four 3-sigma measurements or seven 
2-sigma measurements). We summed the total number 
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of expected events over the 5000 stars and then scaled 
those results to 3.5 years of observation of 150,000 stars, 
assuming that 25% of these stars will be identified as 
variable and not be usefuljl]. That is we assume 390,000 
star-years. Our results are given in Figure 1. In order 
to turn these results into the potential sensitivity of de- 
tecting PBH dark matter, we calculated the 95% C.L. for 
each PBH mass, assuming that no events were detected 
and that there was no background. These potential lim- 
its are shown in Figure 1, along with limits from earlier 
experiments. We found very little difference between the 
two S/N requirements mentioned above, so show only the 
case for 4 sequential 3-sigma event selection. 

We see from Figure 1 that we have the potential to de- 
tect or rule out PBHs as the primary constituent of DM 
over the mass range 5 x lO'-^^Af© < mpBH < 10"'* M©. 
Current hmits from the MACHO/EROS experiments [u 
are also shown as a dashed line. There are also limits 
ruling out halo fractions, 1, from femtolensing of 

gamma ray bursts (GRB)[l9|, but they run from about 
IO-i^Mq < mpBH < 10-i3Mq, off to the left of Fig- 
ure 1, and other limits from picolensing of GRBs ruling 
out / < 4 from 4 x W^^^Mq < mpBH < 8 x 10'^° Mq, 
are in our mass range, but too weak to be seen in our 
plot. We see that a microlensing search for PBHs through 
Kepler data has the potential to extend the mass sensi- 
tivity by almost two orders of magnitude below the MA- 
CHO/EROS limits which exclude DM masses down to 
around 2 x 1O~*M0. Note that commonly quoted (e.g. 

limits from EROS alone 12 1 exclude masses down to 
around 6 x 10~^AfQ and are not as strong as the earlier 
combined MACHO and EROS limits. There are no other 
limits in the mass range just below 2 x 1O~^M0, so the 
capability of Kepler to search for these PBHs is unique. 



DISCUSSION AND FUTURE WORK 
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FIG. 1. Top panel: The expected number of events, A^cxp (de- 
fined as 4 sequential measurements with flux 3-sigma above 
average) in 390,000 star-years of Kepler data. The thin hor- 
izontal line shows A'^exp ~ 3, the limit for a 95% C.L. if no 
events are detected. Bottom panel: Potential 95% C.L. ex- 
clusion of PBH dark matter. The area above the thick line 
would be ruled out if no events were seen in 390,000 star-years. 
Also shown as a thin dashed hue is the hmit from combined 
MACHO/EROS LMC microlensing searches[lH- The thin 
horizontal line shows a halo consisting entirely of PBHs with 
local DM density p = 0.3GeVcm"^. 



In this theoretical paper we suggest a method to de- 
tect or rule out PBH DM over a large and unexplored 
mass range. If nothing is detected the method has the 
potential to rule out around 40% of the total remaining 
mass range for PBH DM. While we emphasized PBHs, 
since microlensing depends only upon lens mass and size, 
such an experiment would detect or rule out any massive 
compact halo object dark matter in the mass range de- 
scribed. 

The next step is clearly to perform the analysis sug- 
gested here, and we have begun this task using pub- 
licly available Kepler data. This experiment requires a 
deep understanding of the lightcurve data including sys- 
tematic effects and possible backgrounds. It requires a 
careful selection of microlensing candidates and an accu- 
rate calculation of the efficiency of that selection method, 
as well as an understanding of false positives and back- 
ground events. Since Kepler transit searches have already 



turned up many stellar flares on M and K dwarfs [20[ us- 
ing selection criteria similar to those used above, a good 
method of distinguishing microlensing from stellar flares 
will be needed. However, there are some powerful dis- 
criminants against backgrounds. The lightcurves should 
match a limb-darkened finite source microlensing shape, 
and if enough events are found the rates should exhibit 
the strong dependence on distance and source radius pre- 
dicted by the above formulas. 

Several effects in the data may mean that the actual 
limits, supposing no PBHs are detected, will be different 
from what we display in Figure 1. First, our simplis- 
tic removal of 25% of the stars as variable Q needs to 
be replaced by detailed analysis of all the lightcurves. 
Our calculation of thresholds using Gaussian statistics 
needs redoing in light of actual lightcurve data and other 
sources of noise that exist in real data. A more careful 
treatment of stellar limb darkening may shorten the aver- 
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age duration of events. The actual selection of candidate 
events may turn-up poor S/N microlensing candidates 
that will reduce the significance of any limits derived. 
These effects will reduce the expected number of events 
used to calculate Figure 1. Of course, if the Kepler mis- 
sion runs longer than 3.5 years the sensitivity would be 
better than that in Figure 1, and it also might be possi- 
ble to use lightcurve shape discrimination to detect lower 
magnification events, thereby increasing sensitivity. 

There are also several assumptions and simplifications 
in our analysis that could change the results. For example 
our estimate of the distances to the source stars is overly 
simple, as is our treatment of limb-darkening. We ignored 
the MW halo model and any transverse velocity of the 
Sun and Kepler source stars. Since the Kepler field is 
in the direction of the Sun's orbit about the MW, any 
rotation of the DM halo in the direction of the disk would 
give support to the halo and not contribute to transverse 
microlensing velocity. This would mean that the average 
transverse velocity of the PBHs relevant for microlensing 
might be lower, implying the average duration of events 
might be longer. Longer events are easier to detect, so if 
this is the case, our limits might be conservative. 
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